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1.  Introduction 


The  technique  of  photoluminescence  (PL)  is  customarily  used  to  investigate  band  structure  and  other 
single-particle  states  of  a  solid  by  monitoring  electronic  transitions  induced  by  incident  high-energy 
photons.^  However,  PL  excitation  also  leaves  the  solid  in  a  high-energy  state  (here  referred  to  as  a 
solid-state  plasma  (SSP)  characterized  by  various  many-body  correlated  subsystems.  Among  these 
are 


•  the  electron  and  hole  gases,  consisting  of  free  electrons  (denoted  by  e)  and  free  holes  (denoted 
by  h) 

•  the  exciton  gas,  consisting  of  bound  electron-hole  pairs  (denoted  by  eh) 

•  the  biexciton  gas,  consisting  of  pseudo-hydrogen  molecules  (denoted  by  eehh) 

•  the  trion  gases,  consisting  of  bound  carrier  triplets  {eeh  denotes  an  e-trion,  ehh  an  /2-trion) 

•  electron-hole  droplets 

Because  all  of  these  species  are  present  at  once  when  a  solid  is  under  PL  excitation,  their  mutual 
interactions  and  carrier  kinetics  are  important  in  interpreting  experimental  data. 

In  3-dimensional  (3-D)  systems  (i.e.,  bulk  solids)  SSPs  are  usually  observed  only  at  low 
temperatures,  since  the  constituent  carriers  consist  of  weakly  bound  “molecules”  made  up  of  the  PL- 
excited  free  electrons  and  holes.  However,  these  entities  have  much  larger  binding  energies  in  2- 
dimensional  (2-D)  systems,  notably  in  the  transition-metal  dichalcogenides  (TMDs)  such  as 
molybdenum  disulfide  (M0S2)  and  tungsten  diselenide.  High-order  correlated  states  (excitons, 
trions,  biexcitons,  triexcitons,  etc.)  have  been  observed  experimentally  in  PL-excited  carrier  gases 
primarily  in  pseudo-2-D  systems  such  as  quantum  wells^;  however,  the  recent  observation  of  trions 
in  M0S2  by  Mak  et  al.^  has  generated  renewed  interest  in  these  truly  2-D  systems,  and  in  particular 
has  motivated  the  writing  of  this  report. 


2.  Formalism 


It  is  instructive  to  examine  SSPs  using  lowest-order  kinetic  theory,  which  is  formulated  in  terms  of 
systems  of  spatially  uniform  rate  equations."^  A  system  of  this  kind  usually  consists  of  N  nonlinear 
differential  equations,  one  for  each  of  the  N  distinct  carrier  species,  with  N  first-order  time 
dependent  densities: 

(t)  =  G. -91. (/p/2,/3,. ..,/^)  ,  i  =  l,2,3,...,A  (1) 
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where  are  source  terms  originating  from  the  incident  PL  photon  flux  and  denotes  the 

recombination  rate  for  the  i  -  th  species.  In  this  context,  “bimolecular”  recombination  implies  that 
the  consist  of  sums  of  quadratic  terms,  each  of  which  is  the  product  of  2  carrier  densities. 

Positive  terms  in  tend  to  increase  the  density  of  carrier  z,  negative  terms  to  decrease  it.  Currently, 

the  most  interesting  PL-induced  systems  are  the  2-D  monolayers,  in  which  recombination  dynamics 
can  be  complicated;  for  simplicity’s  sake  the  modeling  described  in  this  technical  report  is  based  on 
simple  (i.e.,  “pure”)  bimolecular  recombination  laws,  which  can  be  rigorously  justified  for  3-D 
systems. 


Whereas  a  linear  time-dependent  system  of  first-order  differential  equations  has  only  trivial  steady- 
state  solutions  (all  carrier  densities  zero),^  nonlinear  systems  of  this  kind  can  have  multiple  nonzero 
steady-state  solutions  given  by  the  equations 


=  0  , 


z  =  l,2,3,...,iV 


In  this  technical  report,  6  types  of  carriers  are  considered  as  listed  in  Table  1. 


(2) 


Table  1  Densities  of  carriers  in  a  PL-excited  solid-state  plasma 


Carrier  Species 

Density  Symbols 

Electrons 

Ve 

Holes 

Vh 

e-Trions 

Te 

/z-Trions 

Th 

Excitons 

n 

Photons 

_ 2 _ 

3.  Reactions 


The  time  dependence  of  each  carrier  species  is  governed  by  a  set  of  reactions,  which  defines  a  path 
in  the  phase  space  of  the  system  that  leads  to  various  steady  states.  Some  of  these  reactions  convert 
one  member  of  a  species  to  a  different  member  of  the  same  species,  which  can  change  the 
temperature  of  the  carrier  gas  but  not  its  density;  this  process  can  act  as  a  kind  of  “catalysis”  that 
changes  other  carrier  densities  while  leaving  the  catalyst’s  density  unchanged. 

Let  us  examine  each  of  the  reactions  in  Sets  A-E  in  detail: 

A.  Reactions  that  change  the  electron  density. 
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Note  that  Reaction  3  in  Table  2  is  catalytic  as  defined  above,  because  it  changes  the  electron  and 
e-trion  densities  without  changing  the  density  of  excitons.  For  information  on  Reaction  4  in  this 
table  see  Crowne  et  al.^  and  Blythe  and  BloorJ  The  rate  equation  for  electrons  is 


dvp  1  2 

- =  G-  AvpV,  — Px,  Vp  +  YpTiXp  +  Jn 

Qf  ^  h  2  ^ 


(3) 


where  G  is  the  generation  rate  due  to  the  incident  photon  flux  from  the  PL  laser,  which  is  fixed. 


Table  2  Reactions  that  change  Ve 


Reaction 

Rate  Constant 

1. 

Ve  +  Vh  ^  (p 

A 

2. 

Ve  +  Xh  Ex 

P 

3. 

Ex  +  Te  Ve  +  Ex’ 

Ye 

4. 

Ex  +  Ex’  Ve  +  Vh 

J 

B.  Reactions  that  change  the  hole  density. 

Note  that  Reaction  3  in  Table  3  is  also  catalytic,  changing  the  /t-trion  and  hole  densities  without 
changing  the  density  of  excitons.  The  rate  equation  for  holes  is 

dv,  1 

=  G-  AvgVj^  -—NxgV^  +  +  Jn  (4) 

The  factors  of  appearing  in  rate  Eqs.  3  and  4  serve  to  “balance”  these  equations  in  the  chemical 
sense  (i.e.,  each  exciton  produced  requires  2  parent  particles). 


Table  3  Reactions  that  change  Vh 


Reaction 

Rate  Constant 

1. 

Ve  +  V/j  ^  (p 

A 

2. 

V/j  +  Te  Ex 

P 

3. 

Ex  +  T/j  ^  Vh  +  Ex’ 

Yh 

4. 

Ex  +  Ex’  Ve  +  Vh 

J 

C.  Reactions  that  change  the  e-trion  density. 

The  third  reaction  is  listed  in  both  Table  2  and  Table  4  because  it  alters  the  densities  of  both 
e-trions  and  excitons,  as  the  rate  equation  below  shows: 


dxp  1 

—  =  Lnvg  --NxgV^  -  Ygnxg  -  K^x^x^ 


(5) 
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Table  4  Reactions  that  change  Te 


Reaction 

Rate  Constant 

1.  Ex  +  Ve  Xe 

L 

2.  vh  +  Xe  Ex 

N 

3.  Ex  +  Xe  Ve  +  Ex’ 

Ye 

4.  Xh  +  Xe  -Y  Ex 

Ke 

D.  Reactions  that  change  the  /t-trion  density. 


Again,  the  third  reaction  appears  in  both  Table  3  and  Table  5  because  it  alters  the  densities  of  both 
/t-trions  and  excitons,  as  the  rate  equation  below  shows: 


1 

— —-Mnv, Px,Vf,-Y,nx,  -K,Xf,x, 

Qf  n  2  ^  n  n  n  n 


(6) 


Tables  Reactions  that  change  Ta 


Reaction 

Rate  Constant 

1.  Ex  +  v/j  xh 

M 

2.  Ve  +  xh  Ex 

P 

3.  Ex  +  T/j  ^  vh  +  Ex’ 

Yh 

4.  Xe  +  Xh  Ex 

Kh 

E.  Reactions  that  change  the  exciton  density. 

The  second  and  third  reactions  are  catalytic  agents  for  the  loss  of  excitons.  This  implies  that  the 
constants  R  and  S  will  appear  only  in  the  exciton  rate  equation.  The  overall  rate  equation  for  excitons 
is 


dn 

dt 


Q^e^h  ) " 


(V) 


Reactions  that  change  n  are  shown  in  Table  6. 
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Table  6  Reactions  that  change  n 


Reaction 

Rate  Constant 

1. 

Ve  +  Vh  Ex 

e 

2. 

Ex  +  Ve  -N 

Ve’ 

R 

3. 

Ex  +  V/!  -N 

Vh’ 

S 

4. 

Ex  +  V/!  -N 

Xh 

M 

5. 

Ex  +  Ve  -N 

Xe 

L 

6. 

Ve  +  Xh  -N 

Ex 

N 

7. 

Vh  +  Xe  -N 

Ex 

P 

8. 

Ex+Ex’  -Nve  +  vh 

J 

9. 

Xe  +  Xh  -N 

Ex  +  Ex’ 

Ke  +  Kh 

4.  Charge  Neutrality  Ansatz 


An  ansatz  that  is  commonly  used  in  the  theory  of  photoeonduetivity  to  solve  equations  of  this  kind  is 

loeal  eharge  neutrality,  i.e.,  partieles  are  ereated  and  annihilated  in  pairs  so  that  there  is  no 

aeeumulation  of  eharge.  However,  this  approaeh  is  rigorously  eorreet  only  for  a  system  with  2 

earrier  speeies.  An  appealing  way  to  generalize  it  to  our  ease  is  to  assume  it  applies  both  to  eleetrons 

and  trions  separately,  in  whieh  ease  =  v,  =v  and  =  x,  =  x  • 

^  h  ^  h 


Let  us  eonsider  the  full  system  of  5  equations  under  the  strong-neutrality  eondition.  As  a  further 
simplifieation,  let  us  set  the  trion-trion  eonstants  Ke  and  Kh  to  zero.  Then  the  rate-equation  system 
beeomes 


dv  9  1 

—  -  G-  Av^  — Ptv  +  Ypnx  +  J 

dt  2  ^ 

dx  1 

—  =  Mnv  — Pxv  -  Y,  nx 

dt  2  h 

dv  1  \  9 

—  -  G-  Av^  — Nxv  +  Y,  nx  +  Jn 

dt  2  h 

dx  1 

dt  2  ^ 


dn 

dt 


Qv^  -[R  +  S  +  M  +  N)vn 
+  (T’  +  A)tv-  2 


(8) 


Beeause  there  are  5  equations,  this  assumption  over  determines  the  system.  However,  it  is  easily 
shown  that  for  the  speeial  ease  where  N  =  P  =  Il  ,L  =  M  =  A  ,Yq  =  Yy^=P  ,  and  R  -  S  =  Y  ,  these 

equations  take  the  form 
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(9fl)  ^  =  G-Av^--Uxv  +  Tm  +  Jn^ 
dt  2 

dx  1 

(9b)  —  =  Anv — riiv-Twi 
dt  2 

dv  7  I  9  (9 

(9c)  —  =  G-Av^ — UTv  +  Tnx  +  Jn^ 
dt  2 

dx  1 

(9d)  —  =  Anv — riiv  -  T«t 
dt  2 

(9e)  ^  =  Qv^ -2{i:  +  A)vn  +  2Uxv-2Jn^ 

Note  that  Eqs.  9a  and  9c  are  now  identical,  as  are  Eqs.  9b  and  9d,  respectively.  Then  in  this  limit  the 
system  size  is  reduced  from  5  to  3  equations  and  the  number  of  rate  constants  drops  from  10  to  5,  so 
that  for  this  special  set  of  material  constants  the  charge  neutrality  method  yields  a  valid  solution,  at 
least  formally. 

To  find  the  steady-state  densities,  we  set  the  time  derivatives  to  zero,  which  reduces  the  system  (Eq. 
9)  to  the  following  set  of  3  algebraic  equations: 


(lOa)  ^Tlvx-Ynx  =  G-Av^  +Jn^ 

(lOh)  ^nvx  +  Ynx  =  Anv  (10) 

(10c)  2{Y  +  A)nv  +  2Uxv  =  Qv^-2Jn^ 

By  setting  v  =  xn  and  x  =  yn,wQ  obtain  the  following  solutions  to  Eq.  10: 


7 

7 


A  + 


U-2Yx 

n  +  2Yx 


Ax-—Jx^ 

2 


2Ax 


X  r  =  - V  r 

J  U  +  2Yx  J 

where  x  is  a  real  positive  root  of  the  cubic  polynomial  dAflx  =  ijx^  +  2Ax  -  0](n  +  2Yx) . 


(11) 


To  proceed,  we  need  numbers  for  the  5  rate  constants.  To  get  these  constants,  we  can  use  the 
Eangevin  method^  (see  Appendix),  which  partitions  the  recombination  into  a  “hunting”  process 
where  the  reactants  approach  each  other  in  real  space,  followed  by  an  “on-site”  event  of  short 
duration.  Eangevin  argued  that  since  the  contribution  of  the  latter  process  is  small,  most  of  the 
recombination  time  comes  from  the  “hunting”  phase.  We  will  argue  that  that  this  scenario  is 
appropriate  for  recombination  kinetics  in  TMDs  for  the  following  reasons: 
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1)  The  large  binding  energy  of  excitons  and  other  many-body  complexes  in  TMDs  suggests  that 
they  are  compact  in  real  space  (i.e.,  Frenkel-like  objects)  with  large  effective  masses  and  hence 
low  mobilities,  a  characteristic  that  is  enhanced  by  strong  electron-phonon  coupling,  poor 
material  quality,  etc. 

2)  The  J-state  nature  of  the  conduction  and  valence  bands  of  TMDs  implies  that  for  such 
complexes  recombination  takes  place  between  atomic  states.  This  implies  that  the  matrix 
elements  for  atomic-level  decay  processes  are  large,  and  hence  the  “on-site”  recombination 
time  is  small. 

Interaction  between  pairs  of  complexes  is  governed  by  the  electrostatic  forces  between  them. 
Although  the  interaction  mechanism  for  recombination  is  strictly  quadratic  in  this  technical  report. 
Table  7  lists  the  various  force  laws  that  mediate  these  interactions  and  enter  into  Langevin’s  theory. 
We  will  treat  excitons  as  neutral  polarizable  entities  whose  interactions  with  charge  carriers  and 
trions  resemble  those  of  nonpolar  molecules. 

Table  7  Force  laws  that  mediate  Langevin  recombination  between  carrier  species:  |Ue,/!— electron/hole  mobilities, 
|Li,^_trion  mobilities,  and  |Uex— exciton  mobilities. 


Recombining  Pair 

Rate  Constant 

Force  Law 

Langevin  Coefficient 

1.  Ex,  V 

Nonpolar  dipole  +  charge 

A,  2 

2.  V,  T 

Charge  +  charge 

n 

r-2 

3.  Ex,  T 

Nonpolar  dipole  +  charge 

Y 

~~[^ex  V-tr) 

SSq 

4.  T,  T 

Charge  +  charge 

K 

r-2 

SSq 

5.  Ex,  Ex 

Nonpolar  dipole  +  Nonpolar  dipole 

J 

r-7 

London  dispersive 
force 

5.  Experimental  Orders  of  Magnitude  for  Rate  Constants 


In  order  to  test  the  relevance  of  these  calculations,  we  can  estimate  the  densities  of  the  various 
carrier  species  predicted  by  this  theory  within  a  2-D,  PL-excited  M0S2  monolayer.  In  our  steady- 
state  experiments,^  we  used  a  532-nm  laser  (photon  energy  2.33  eV)  to  create  an  incident  intensity  of 
140  pW/pm^  at  the  sample  surface,  which  corresponds  to  a  photon  flux  of 
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~3.75  X 10^^  photons/ *  s j .  A  quantum  yield  for  electron-hole  pair  creation  of  0.004  would 
create  a  generation  flux  G  within  the  monolayer  of  '-'1.5  x  10^®  eh  pairs/ ^cm^  *  sj .  Assuming  that  the 

Langevin  approach  is  appropriate,  we  can  use  our  in-house  samples  of  2-D  M0S2  with  field-effect 
transistor  mobilities  p,  =  10  cm^  /  ( V-s)  to  derive  a -\\.ql  |sSq  j  =  1.8xl0  ^  cm^  /  s  and 

_ A  O 

A  =  Ana  =  2.3  x  10  cm  /s  .  If  only  the  ^f-type  recombination  process  operated,  this  value  in 
tandem  with  G  would  predict  an  equivalent  steady-state  carrier  density 

=  8.07x10^'  cm“^ 

6.  Numerical  Simulation 


By  restoring  the  time  derivatives  to  the  truncated  version  of  Eq.  9,  we  obtain  a  system  of  3  ordinary 
differential  equations  that  is  amenable  to  solution  by  numerical  simulation,  notably  by  using  the 
routine  NDSolve  in  the  Mathematica  package.  Let  us  assume  that  all  of  the  rate  constants  are  the 
same  order  of  magnitude,  as  listed  in  Table  8. 


Table  8  Rate  constants  for  Eq.  2 


Rate  Constant 

Valne 

n 

4.xl0“‘^ 

A 

4.x 10“^ 

Y 

1.x 10“^ 

J 

ExlO”"^ 

Q 

l.xl0“^ 

Restoring  the  time  derivatives  results  in  the  time-dependent  system 

dv  7  1  1  7 

(12a)  —  =  G-Av^ — nxv  +  Ym  +  -Jn^ 

dt  2  2 

1 

(12h)  —  -  Anv — rixv-Ynx  (12) 

dt  2 

(12c)  —  =  Qv^-  lAvn  +  2nxv  -  Jn^ 

dt 

with  the  same  fixed-point  solutions  as  Eq.  1 1 ; 
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-1 


7 

V 

V 


A  + 


n-2Yx 
n  +  2Yx 


Ax  — Jx 
2 


2 


-XVf 

2Ax 

= - ^  f 

n  +  2Yx  / 


(13) 


where  x  is  again  a  real  positive  root  of  cubic  polynomial  4Anx  =  |  Jx^  +  2Ax  -  0j(n  +  2Yx) .  In 

principle  as  many  as  3  steady  states  are  therefore  possible.  Using  the  experimental  numbers 
G  =  1.5xl0^®  ehpairs/|cm^  "sj ,  A  =  2.5  x  10”^ ,  and  ^4  -  2.3x  10~^cm^ /s  together  with  the  rate 

constants  listed  in  Table  8,  we  find  that  there  is  one  real  root  for  the  cubic,  which  gives 

V j.  =8.16x10^^  cm~^ 

=8.49x10^*’  cm"^  (14) 

Ty  =1.69x10^  cm~^ 

For  these  values  of  the  various  rate  constants,  in  steady  state  the  majority  of  the  carriers  are  eh  pairs, 
followed  by  excitons  and  then  by  trions.  Running  the  time  simulation  by  switching  on  the  laser  at  t  = 
0  with  starting  values  v{t  =  O)  =  =  O)  =  T(t  =  O)  =  0  gives  the  plot  shown  in  Fig.  1; 


LogN 


Fig.  1  Time  dependence  of  carrier  species  for  A  =  4.  x  10  Red — electron-hole 

pairs,  magenta — excitons,  and  blue — trions. 

Note  that  the  fastest  rise  time  is  for  the  eh  pairs,  followed  by  excitons  and  then  by  trions. 
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It  is  also  noteworthy  that  the  constant  A,  which  converts  an  exciton  into  a  trion,  is  very  small. 
Suppose  we  increase  it  to  2.5  x  1 0  ^ .  The  new  simulation  results  shown  in  Fig.  2  exhibit  several  new 
features.  First  of  all,  the  fixed  point  has  moved; 

V j.  =8.16x10^^  cm~^ 

=3.99x10^2  ^m“^  (15) 

Ty  =6.41x10^^  cm~^ 

and  now  predicts  that  there  are  now  more  trions  than  excitons.  In  addition,  the  simulation  now  shows 
changes  in  the  transient  behavior: 

1)  The  eh  density  seems  to  saturate  around  10  ns,  but  then  develops  a  minimum  in  the  range  from 
10  ns  to  140  ns  before  reaching  a  true  steady  state. 

2)  Initially  there  are  more  excitons  than  trions,  with  a  crossover  around  20  ns.  The  exciton  curve 
is  nearly  linear  over  a  wide  range  from  10  ns  to  50  ns. 


Fig.  2  Time  dependence  of  carrier  species  for  A  =  4.  x  10^.  Red — electron-hole 
pairs,  magenta — excitons,  and  blue — trions. 


These  features  are  more  pronounced  when  A  =  2.5  x  10 


-3 


,  for  which  the  new  fixed  point  is  at 
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(16) 


V j.  =  8.16 X 10^^  cm 
tij.  =1.15x10^^  cm“^ 

Xy  =1.35x10^'^  cm~^ 

Next,  Fig.  3  graphically  exhibits  the  time  dependence  of  carrier  species  for  A  =  4.  x  10^^ 


Fig.  3  Time  dependence  of  carrier  species  for  A  =  4.  x  10  .  Red — electron-hole 
pairs,  magenta — excitons,  and  blue — trions. 

Another  parameter  of  interest  is  J,  which  controls  the  decay  of  excitons  into  eh  pairs  by  exciton- 
exciton  collisions.  This  reaction  is  complementary  to  the  parameter  ^4,  which  mediates  the  ereation 
of  excitons  from  photons.  For  the  simulation  plotted  inFig.  l,Jwas  set  to  l.Ox  10”^;  deereasingj 
to  1 .0  X  10”^  results  in  a  new  fixed  point  and  the  transient  behavior  shown  in  Fig.  4: 

V J.  =8.16x10^^  cm~^ 

Uy  =4.89x10^^  em“^  (17) 

Xy  =7.52x10^®  cm“^ 
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LogN 


Fig.  4  Time  dependence  of  carrier  species  for  J  =  4.  x  10  Red — electron-hole 
pairs,  magenta — excitons,  and  blue — trions. 

Decreasing  /  results  in  a  rapid  rise  time  and  saturation  for  the  eh  pairs,  while  excitons  and  trions 
increase  more  slowly,  with  the  former  saturating  above  the  eh  pairs  and  the  latter  below  them. 

Increasing /to  1.0  x  10“^  results  in  the  fixed  point 

V j.  =  8.16x10^^  cm~^ 

=8.49x10^°  cm“^  (18) 

Ty  =1.69x10^  cm~^ 

and  the  time  history  shown  in  Fig.  5.  Now  all  of  the  carrier  densities  rise  and  saturate  rapidly,  with 
the  trion  population  the  smallest  and  eh  pair  population  the  largest. 
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LogN 


Fig.  5  Time  dependence  of  carrier  species  for  J  =  4.  x  10  .  Red — electron-hole  pairs, 
magenta — excitons,  and  blue — trions. 


The  final  parameter  of  interest  is  Q,  which  controls  the  creation  of  excitons  by  electron-hole 

collisions.  It  is  noteworthy  that  all  the  simulation  results  presented  up  to  now  have  predicted  the 

same  steady-state  eh  pair  density,  i.e.,  8.16  xlO"  cm~^  ,  regardless  of  how  the  other  reaction 

constants  are  varied.  To  investigate  this  behavior,  let  us  use  the  system  (Eq.  8)  to  construct  the 

quantity  2  — +  2—  +  —  -  This  leads  to  a  detailed-balance  type  of  equation  of  the  form 
dt  dt  Bt 


2 


dt 


+  2  —  +  —  =  2G-i2A-Q)v^ 
dt  dt  ^  ’ 


(19) 


If  we  postulate  that  the  system  has  a  steady  state,  we  can  state  that  the  combination  of  derivatives 

Btv  B  B  I  2G 

2  —  -1-2  —  -r  —  =  o;  then  (19)  implies  the  relation  v  =  I - .  For  the  constants  we  have  chosen, 

Bt  Bt  Bt  ]j  2A  —  Q 

this  gives  8.16  xlO"  cm~^  ,  i.e.,  in  exact  agreement  with  the  value  of  v  predicted  by  the  simulations. 
This  result  would  be  useful  except  that  the  argument  is  actually  circular,  in  the  sense  that  if  we 
consider  the  case  where  Q>2A,  the  system  clearly  does  not  have  a  steady  state.  In  fact,  as  shown  in 
Fig.  6,  the  response  predicted  by  simulation  is  “explosive”,  and  the  simulation  program  responds  to 
this  state  of  affairs  by  terminating  at  a  time  of  approximately  60  ns,  beyond  which  it  cannot  solve  the 
problem. 


13 


Fig.  6  Explosive  behavior  of  simulation  with  Q  >  2A.  System  (Eq.  8)  has  no  steady 

state,  forcing  the  simulation  program  to  terminate  at  a  time  of  approximately  60  ns. 


7.  Strong  vs.  Weak  Neutrality 


On  the  surface,  the  assumption  of  equal  densities  should  be  a  reasonable  starting  point  for  analysis  of 
the  carrier  dynamics.  However,  since  the  correct  problem  involves  5  rate  equations,  there  are  2  ways 
to  implement  charge  neutrality:  the  ansatz  of  Section  4  above,  which  we  will  refer  to  as  “strong 
charge  neutrality”,  and  also  a  “weak  charge  neutrality”  condition,  for  which  Vg  +  Tg  =  .  This 

approximation  reduces  the  number  of  equations  from  5  to  4.  Which  is  the  right  choice? 


One  way  to  test  the  validity  of  the  strong-neutrality  condition  is  to  check  for  stability  of  the  system 
when  this  condition  is  weakly  violated.  For  this  we  use  the  full  5 -equation  system  with  the  special 
set  of  coefficients  but  allowing  all  5 -carrier  densities  to  be  unequal,  which  reads 


dvp  1  9 

=  G-  4vgV^  -  — Hx^Vg  +  Y«Tg  +  Jn 


dt 


1  7 

=  G  -  AvpV, - Hx^v,  +Y«x,  +  Jn 

^  n  2  ^  " 


dt 

dxp 


dt 

dx 

~dt 

dn 


=  A«Vg  -  — ffXgV^  -  Y«x^ 


^  =A«v^-^nx^Vg-Y«x^ 


^  =  Qv,v^  - 2A(vg  +  vj «  +  n (xgv^  +^h^e)- 


(20) 
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The  steady-state  solution  to  this  system  satisfies  the  strong-neutrality  eondition;  let  us  perturb  the 
densities  around  this  solution  as  follows: 


Vg  =  V  +  6vg 

V,  =v  +  5v, 
h  h 

Tg  =  T  +  SXg 

T,  =t  +  5t, 
h  h 

«  =  «Q  +  6« 


This  leads  to  the  linear  perturbation  equation 


'5vg^ 

^5vg^ 

5xg 

=  M 

5xg 

V  y 

with 


M  = 


(  1  ^ 

fiv  +  — ITt 

I  2  ) 

-{Av) 


-(fiv) 


{  1  ^ 

fl  ) 

■  Av+-nx 

-Uv 

1  2  ; 

U  ) 

-nv 

v2 


Tx  +  2Jn, 


0 


Yx  +  2Jn, 


0 


‘0 


-nx 

u 


f  1 

-nv  +  Yn, 
[2  ' 


fl  ^ 

fl  f 

-Ux 

U  ) 

A«o 

0 

-Uv  +  Ym. 

U  Oj 

Qv- 

V 

2a«q  +nx 

Qv  —  2A«q  +  fix 

Uv 

Uv 

(21) 


(22) 


(23) 


Av  -  Yx 

Av- Yx 
+  Avj 

The  stability  eriterion  is  that  the  eigenvalues  of  this  equation  all  have  negative  real  parts;  this  ensures 
that  any  departure  from  the  fixed  point  will  deeay  baek  to  it.  Using  the  first  set  of  rate  eonstants  and 
the  eorresponding  fixed  point  of  the  equations 

-3 


V j  =  8.16x10^^  em 


n j  =8.49x10^*^  em  ^ 
Ty  =1.69x10^  em~^ 


(24) 


the  eigenvalues  A,/  are  easily  evaluated  numerieally: 
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=-3.81x10^ 

^2  --1.79x10^ 

^3  =-1.68x10^  (25) 

^4  --1.13x10^ 

^3  =+1.15x10'^ 

The  sign  of  reveals  that  in  faet  the  solution  to  the  system  with  strong  neutrality  is  not  a  fixed 

point  of  the  system,  from  whieh  we  may  eonelude  that  the  true  system  has  a  steady  state  with  weak 
neutrality  (i.e.,  the  SSP  has  more  eleetrons  than  holes  and  more  A-trions  than  e-trions  to  eompensate 
it). 

8.  Conclusions 


The  analysis  given  here  only  seratehes  the  surfaee  of  this  problem,  both  mathematieally  and  with 
regard  to  the  physies  of  transport  in  2-D  systems.  Use  of  Langevin’s  formula  should  elearly  be  a 
point  of  eontention;  the  work  of  Juska  et  al.^*^  on  low-mobility  layered  systems  for  solar  eells 
indieate  that  the  dimensionality  of  the  system  enters  strongly,  eonverting  the  simple  bimoleeular  law 
in  3-D  to  a  2-D  expression  involving  fraetional  powers  of  the  densities.  In  addition,  Greenham  and 
Bobbert^^  have  observed  that  the  transport  is  partially  diffusive;  that  is,  the  Langevin  dependenee 
takes  into  aeeount  only  the  drift  of  the  partieles  in  a  Coulomb  field.  Finally,  it  is  possible  to  take 
some  temperature  dependenee  into  aeeount  by  introdueing  the  “Coulomb  radius” 

=  q  I  AmE^kT  as  a  boundary  between  Langevin’s  earrier  approaeh  and  the  “blaek  hole”  where 

the  earriers  reeombine.  We  plan  to  diseuss  these  matters  in  a  future  publieation. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


M0S2 

molybdenum  disulfide  (a  TMD) 

PL 

photolumineseenee 

SSP 

solid-state  plasma 

TMD 

transition-metal  diehaleogenide 

2-D 

2-dimensional 

3-D 

3 -dimensional 
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